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ABSTRACT 

In the present study, a rigorous heat transfer model 
has been developed for a rotary kiln used for drying and 
preheating of wet solids. The results have been checked 
by applying this model to the non-reacting zone of a cement 
rotary kiln mentioned in a previous work . 

A rotary kiln consists of a refractory lined cylin- 
drical shell mounted at a slight incline from the horizontal 
plane. The kiln is very slowly rotated about its longitudinal 
axis. The wet solid is fed into the upper end of the cylinder 
and during the process is dried and preheated by the counter- 
current flow of the hot gas. Finally, it is transferred to 
the lower end where it is discharged. 

In the present work, a heat transfer model is developed 
which includes radiation from the hot gas and transient 
conduction in the wall and mass and energy balances of the 
hot gas and the solid. The main objectives of this investi- 
gation are to determine the length of the kiln required to 
dry and preheat the wet solids to a desired temperature and 
also the axial gas and solid temperature distribution. 

To compute the radiation heat transfer, the wall 
and the solid surface are divided into twenty elements that 



xi) 

form a gray enclosure. Hot gas which emits, absorbs and 
transmits radiation fills the enclosure. The shape factors 
between the elements are determined and the net radiation 
exchanged is computed using gas-radiation theory. Heat 
conduction in the refractory wall is determined using 
the finite difference technique and assuming heat conduction 
only in the radial direction. Convection heat transport 
due to kiln rotation is included in the energy balance of 
th^ wall elements . The rate of vaporisation of the wet 
solids is determined by making a mass and energy balance of 
the^ solids . 

Although the results show that the predicted length 
is much lower than that of the actual kiln, the axial 
solid and gas temperature distribution Vs. percent kiln 
length matches well with the actual kiln profile. It is 
also seen that there is virtually no circumferential variation 
of refractory wall temperature. The radial wall temperature 
variation is seen to be almost linear. 
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Chapter 1 


1 .1 Introduction 

Although rotary kilns are used in many chemical 
and metallurgical industries, the information that have 
been published regarding the heat transfer analysis and 
design of these units are very few. In the present 
study, a heat transfer model is developed for the non- 
reacting zone of a cement rotary kiln. Basically, 
this model can be used for any rotary kiln which is 
used for drying and preheating of wet solid, e.g., 
preheating of wet iron ore in metallurgical industry. 

A rotary kiln consists of a refractory lined 
cylindrical shell mounted at a slight incline from the 
horizontal plane. The kiln is very slowly rotated 
about its longitudinal axis. The solid is fed into 
the upper end of the cylinder and during the process 
is transferred towards the lower end where it is dis- 
charged. Hot gas is fed into the lower end of the 
cylinder so that the flow is counter current. Fig. 1 
shows the schematic diagram of a rotary kiln. 

In the present study, the kiln can be considered 
to consist of three sections. In the first: section, 
the wet solids are heated to the boiling point of the 
entrained liquid. In the second region, the liquid 
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FIG-1 SCHEMATIC DIAGRAM OF A ROTARY 
KILN 
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evaporates at constant temperature until the feed is 
completely dry. In the third region, the solids are 
heated to some desired temperature and then are dis- 
charged from the kiln. The influence of simultaneous mass 
transfer and 

/chemical reaction has not been considered in this study. 

The efficiency of heating in rotary kiln heaters 
is better because of improved heat transfer to the 
solids. The rotation of the kiln continuously mixes 
the solids and therefore, the heat transfer rate is not 
limited by solid conduction. In addition, the heat is 
transferred from the refractory wall to the solid in 
the region of contact. Due to rotation, the refractory 
wall is alternately heated and cooled by the hot gas 
and the solid material. Fig. 2(a) shows the nomenclature 
of the kiln. The angle, a, is called the fill angle which 
shows the region of volume containing the solid. The 
angle, ©,is measured from the origin in clockwise 
direction. The origin is represented by the line joining 
the center of the kiln and the right extreme point of 
the solid cross section. The internal radius of the 
kiln is represented by the symbol, R. 

Figure 2 (b) shows the heat transfer mechanisms 
occurring in the kiln. The inner wall which is not 
covered by the kiln, the solid surface, and the hot gas 
exchange heat between each other by. radiation. The 
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FIG-2(b) SCHEMATIC DIAGRAM SHOWING HEAT 
TRANSFER PROCESS IN ROTARY KILN 
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solid receives heat, from the inner wall covered by it. 
Also, heat is lost from the kiln outer wall to the 
surroundings . 

1 .2 Literature Survey 

Sass [l] investigated cement kilns and iron ore 
kilns and Manitius, Kureyusz and Kawecki [2] studied 
rotary kilns used for the calcination of basic ammonium 
aluminium sulphate. Sass's work did not consider chemical 
reactions and mass transfer from solid material but 
Manitius et.al. included these effects. Mauson, L. and 
Unger, S. [ 3^ discussed the heat transfer in a rotary 
kiln incinerator and presented the applicable relations. 
Recently Ghoshdastidar et .al [4 ] developed a detailed 
heat transfer model for the rotary kiln waste incinerator. 
This model included radiation between the flame and 
refractory surface, transient conduction in the wall and 
chemical reaction in the waste material. It can predict 
the length of the kiln required to completely bum 
plexiglas (Polymethyl-methacrylate) for various fill 
angles . 

1 .3 Objectives 

The objectives of iihe present work are: 

.(i) to develop a rigorous heat transfer model 
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for rotary kilns used to dry and preheat wet solids, 
e.g. in cement and metallurgical industries. 

(ii) to apply the present model to the rotary 
cement kiln of Sass [l ] and compare the results, i.e. 
predicted length of the kiln and axial solid and gas 
temperature variations. 

(iii) to check the validity of Sass's assumption 
that no circumferential temperature variation exists 
in the kiln. 

(iv) to plot the radial temperature distribution 


in the wall . 
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Chapter 2 

Problem Formulation 


2 . 1 Radiation Heat Transfer between the Hot Gas, the 
Refractory Wall and the Solid Surface 

Heat is exchanged between the hot gas, the inner 
wall of the kiln and the solid surface by radiation. 

Since the temperature of the wall varies circumferen- 
tially and axially, the wall is divided into surface 
elements as shown in Fig. 3. The elements that are not 
covered by the solid and the elements on the surface 
of the solid are exposed to radiation from the hot 
gas. The portion on the wall is divided into 15 
elements and the solid surface is divided into 5 surface 

elements. Convection from the gas is neglected as radiation 
is predominant mode since temperature of the gas is very high 
In this study, for the purpose of calculating 

emissivity and other characteristics, the composition 
has 

of hot gasZbeen taken from Ing. Paul Weber £53, the 
reason being non-availability of gas composition data 
for Sass's kiln. The radiation heat transfer among 
the surface elements of adjoining axial segments is 
found to be insignificant for the segment length of 
importance in this study. Therefore, the radiation heat 
transfer is assumed to occur only among the 20 surface 



Wall surface elements 
Solid surface elements 


FIG- 3 SECTION OF THE KILN SHOWING SURFACE 
ELEMENTS OF THE WALL AND THE SOLID 
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elements in each axial segment. 

The radiation heat transfer is calculated by using 

the theory [6] of gray enclosure containing a gas which emits, 

absorbs and transmits radiation. Thus, if the gas volume is 
black 

enclosed by N^/surfaces at different temperatures, the net 
energy gain for a particular surface can be expressed as 


q . = G , A . - E, . A . 

j 3 b J J 


( 2 . 1 ) 


where G. and E, , are irradiation and black body emission 
j DJ 

per unit area of surface j . Irradiation of surface 
j = Irradiation from gas + Irradiation from all other 
surfaces through gas . Or, 


G.A, = A F .e (T_)E, + Z A. F. . m (T . ) E, . 


j 2 g gj g g bg 


i ij g i bi 


( 2 . 2 ) 


where F^. is the shape factor, e is the emissivity, t 

4 4 

is the transmissivity and and E^ are a T^ and oT^ 
respectively. 

Heat exchange in gray enclosure can be calculated as 


^gray ^black * 2 


for e >0.8 


where e is the emissivity of wall [6l . In our problem emi- 
ssivities of refractory and solid are taken as 0.9 and 0.95 


respectively . 

For the above analysis, shape factors are calculated 
as shown in Appendix A. Gas emissivity at different gas 
temperatures are calculated as given in Appendix B. Solving 
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the above equations (2.1) and (2.2), heat gain/loss of 
each surface is calculated. The temperature of solid 
surface elements are taken to be uniform in each axial 
segment. The refractory wall surface element temper- 
ature is determined by radiation input to the surface 
qj , heat conduction in the wall and the rotational 
speed of the kiln using the method outlined in the 
next section. 

2.2 Conduction Heat Transfer in Refractory Wall 

Since wall elements of the kiln are alternately 

heated and cooled during each revolution by the hot gas 

and the solid respectively, non steady heat conduction 

is present. Only heat conduction in the radial direction 

of the wall is taken into account, assuming negligible 

conduction in the circumferential and axial direction,as the 
circumference and lencrth of the kiln are much greater than the 
thickness of the kiln. For a coordinate system that is fixed to 
rotary kiln, the temperature varies with time and 

consequently energy storage term is needed in the energy 

equation. However, for a coordinate system that is 

stationary, the temperature remains constant and hence 

the heat storage term is not required. However, the 

kiln wall now has a velocity with respect to the coordinate 

system and hence the convective term has to be introduced. 

Because the nonrotating coordinate system is more conveniently 

combined with the radiation boundary condition discussed in 
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the previous section, this coordinate system is chosen 
for the analysis . 

The coordinate system used to model the wall heat 
conduction is shown in Fig. 4. Because the refractory 
thickness is small compared to the diameter, cartesian 
coordinates are used instead of cylindrical coordinates. 
The X and Y coordinates are in the radial and circum- 
ferential directions, respectively. The energy equation 
for the wall then becomes 



rf 



(2.3) 


where U is the circumferential velocity of the wall and 
the thermal diffusivity of the refractory. The 
left hand term of the above equation accounts for the 
energy convection by the refractory wall and the right 
hand term is the radial conduction. 

The boundary condition for the inner surface of 
the refractory is dependent upon whether the surface 
element is exposed to the hot gas or covered with the 
solid. Elements exposed to gas are heated by radiation 
and the heat input is q j , obtained from equation (2.1) 
in the previous section. Convection heat transfer by the 
gas is assumed small compared to radiation and therefore 
is neglected. Surface elements covered by solid exchange 
heat with solid primarily by surface contact. Helmrich and 




FIG-4 SECTION OF THE KILN SHOWING THE 
COORDINATE SYSTEM 
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Schugerl [7] investigated the heat transfer between hi In 
and the solid over a range of particle diameters. They 
expressed heat transfer rates as a product of heat 
transfer coefficient and the difference between the wall 
and bulk particle temperatures. The same approach is 
used in the present study. Therefore, 

= h cs- A J < T sj - T s> (2 - 4 > 

where q. is the heat transfer rate, h is the heat transfer 
J cs 

coefficient, A. is the surface element area, T . is the 
J s J 

refractory wall surface temperature and T g is the solid 
temperature. Since the heat transfer coefficient, h is 

Co 

not available for the non-reacting part of the cement kiln, 

2 

in the present study a value of 64 W/m -K has been used. 

This value is for the calcination range of the kiln as 
given by Helmrich and Schugerl [7 ] . 

The boundary condition for the outer surface of 
the refractory wall is the convection and radiation heat 
losses. N.V. Suryanarayana et al [8] have investigated 
the convection and radiation losses from a rotary kiln 
for various wind velocity, orientation of the kiln, 
rotational speed of the kiln etc. In the present investi- 
gation, convection heat transfer coefficient, h c is 
taken as 10 W/m -K considering wind is blowing over the 
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kiln. Radiation heat transfer coefficient is calculated 
as follows ; 


h = 
r 


a e rf (T - T 1 ) 


(T_ - T* ) 


n 


i.e. h = a £ ^(T + T ) (T^ + T^) 

r rt n « n 00 


(2.5) 


where T r and are outer wall temperature and ambient 

temperature respectively. 

Then combined heat loss from outer surface is 
calculated as 


q j - < h c + V a j (T j - T J 


( 2 . 6 ) 


where T^ is the temperature of outer wall element and 
Aj is the area of that element. 


Equation (2.3) is solved numerically by first 
transforming it into finite difference form. The details 
of finite difference scheme and its method of solution 
are given in Appendix C. 


In the finite difference scheme, the circumferential 
elements of the refractory wall are alligned with the wall 
surface elements as shown in Fig. 3. Five additional 
elements are taken for the region covered with solid. 
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Conduction, through the wall is simulated by sectioning 
the wall into ten layers separated by equal gap. 


2 . 3 Energy Transport in Solid 

As has been mentioned earlier, the kiln can be 
considered to consist of three sections, viz., first 
section in which wet solids are heated to the boiling 
point of the entrained liquid, second section in which 
liquid evaporates at constant temperature and the 
third section in which dry solids are heated to the 
required temperature. 


In the first and third sections, the physical 
processes are similar, the only difference being, in 
the first section it is wet solid and in the third 
section it is dry solid. Fig. 5 shows the energy balance 
on an element of solid contained in an axial segment 
either in the first or in the third section of the kiln. 

From Fig. 5, it can be seen that 

m S.Z C ps T SZ + q l + q 2 = ^S.Z+AZ C ps T S.Z +AZ 

But since no evaporation takes place, we see 

* * 

that m g z = m g z+az - Then above equation can be 
written as follows: 


T = t + 

S.Z+AZ S.Z 


(qi + q 9 ) 


m _ c „ 
S.Z ps 


(2.7) 
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where and q ^ are the total radiation heat transfer 
from the gas and heat transfer from the inner refractory 
wall to the solid respectively and C and m are 

pS O . Zi 

specific heat and mass flow rate respectively of wet 
or dry solid as the case may be. 

In the second section of the kiln, the evaporation 
of the entrained liquid takes place and so the solid 
temperature remains constant at the boiling point of 
the entrained liquid. Fig. 6 shows the mass balance of 
solid in the second section. The mass balance equation 
can be written as 


m 

S.Z+AZ 


m 


S.Z 


"V 


( 2 . 8 ) 


where m^. is the rate of evaporation in that segment. 
Fig. 7 shows the energy balance of wet solid in the 
second section. The energy balance equation can be 
written as 


Vs C ps T S.2 + + «2 = "V (L + C pl T S-Z> + 


+ '"s.Z+iZ C ps T s.z 


Solving equations (2.8) and (2.9), we can see 
that rate of evaporation m^. is obtained as 


(2.9) 
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FfG-5 ENERGY BALANCE OF SOLID ELEMENT 
CONTAINED IN AN AXIAL SEGMENT IN 
THE 1 st OR 3 rd SECTION OF THE KILN 



FIG-6 MASS BALANCE OF SOLID IN THE SECOND 
SECTION OF THE KILN 
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* + q 2 

% = — (2.10) 

where L is the latent heat of evaporation of the entrained 
liquid. 

Since the first 5 elements are on the surface of 
the solid, total radiation heat transfer, q^, is 

5 

= I q 1 (2.11) 

J=1 J 

where q^ is the same as given in equation (2.1). 

Similarly, the heat transfer from the refractory, q 2 , 
is given by 

q 2 = h c« VVi - v - (2 - 12) 


2.4 Energy Balance in Gas 

In the first and third sections of the kiln, only 
heating of wall surface and solid is taking place. Also, 
it is to be noted that the flow of gas is in the opposite 
direction of solid flow. Fig. 8 shows the energy balance 
of hot gas in the first and third sections . 

We can write the energy balance equation as 
^g.Z+^Z C pg T g.Z+ A Z “ m g.z C pg T g.Z + q r (2.13 



y ( L + Cp L 1z) 


^1 



FIG' 7 ENERGY BALANCE OF SOLID IN THE SECOND 
SECTION OF THE KILN 



FIG-8 ENERGY BALANCE OF HOT GAS IN THE 
1 st AND 3 rd SECTIONS 
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where q r is the total heat lost from the hot gas due to 
radiation. But the heat lost by the gas is gained by 
wall surface and solid surface. So 


20 

<3_ = I q, (2.14) 

J=1 j 


where qj is same as given in equation (2.1). 

* • 

is no change in mass flow, m^ z = m^ z+iz * 
equation (2.13), it can be written as 


Since there 
Thus from 


m = t + 

g.Z+AZ g.Z 


Vz C 


pg 


(2.15) 


In the second region, due to evaporation, rate of 
flow is also changing. Fig. 9 shows the mass balance of 
gas in the second section. Referring this figure it can 
be written as 


ra 


g.Z+AZ 


= m 


g.Z 


- \ 


(2.16) 


Again, in the second section, superheating of 
vapour also takes place over and above the heating of 
wall surface and solid. Fig. 10 shows the energy balance 
of gas in the second section. If C^. is the specific heat 
of vapour/ it can be written as 




FIG-9 MASS BALANCE OF GAS IN THE 
SECOND SECTION 



FIG-10 ENERGY BALANCE OF THE GAS IN THE 
SECOND SECTION 
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m g.Z+AZ C pg T g .Z+AZ m - " C — T 


g.z C pg T g.Z + q r + C V 


(T g.Z- T sat 5 


or 


T g.Z+AZ 


m 


g . Z 


<m g.Z - V 


q r + T sat^ 


T „ + — — ^ 


g.Z 


(m g.Z - V C 


pg 


( 2 . 17 ) 


Thus the temperature of gas in the next segment 
can be calculated. 

It should be noted that the temperature of the gas 
and the solid in each axial segment have been taken as 
uniform in this analysis. 
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Chapter 3 
Method of Solution 

A computer program is developed to obtain a numer- 
ical solution to the heat transfer equations. For 
solving the heat transfer problem, the kiln is divided 
into segments of equal length, AZ, in the axial direction. 

In this problem AZ is taken as 1.0 m. Also, the fill 
angle, a, is taken to be 90° (which is same as that of 
Sass's kiln). The solid temperature and mass flow rate 
are known at the inlet of the kiln. The outlet gas 
temperature at the inlet of the kiln, since the flow is 
countercurrent, is also assumed to be known and also 
the rate of flow. The percentage of liquid in the incoming 
wet solid is also known. 

Thus the heat transfer equations can be solved at 
the first segment in the preheating zone. The temperature 
of solid leaving the first segment is calculated using 
equation (2.7). The temperature of gas in the second 
segment is also calculated using equation (2.15). Now the 
second segment can be analysed and it is proceeded 
similarly. 

Once the wet solid temperature reaches the boiling 
temperature of the entrained liquid, evaporation starts. 
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In this section, the mass flow rate of both solid and 
gas and the temperature of gas in the next segment 
are calculated using equations (2.8), (2.16) and 

(2.17), respectively. The temperature of solid remains 
constant in this section. 

Once the whole quantity of liquid is evapora- 
ted, the post heating starts. Here also similar to 
first section, equations (2.7) and (2.15) are used for 
calculating the temperatures of solid and gas respect- 
ively, in the next segments. 

Once the temperature of solid reaches the required 
value, the computation stops and hence the length of the 
kiln required to dry the wet solid and preheat to a 
desired temperature can be calculated. 

Figure 11 shows a concise flow chart of the computer 
program developed to solve the problem. The parameters 
and their values that have been given as input to the 
program are listed in Table 1. A copy of the computer 
program is listed in Appendix D. 




FIG-11 FLOW CHART OF THE PROGRAM 
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Table 1. Input Data to Program 


I . Rotary Kiln 


A. Diameter 

* 

• 

1.95 m 

B. Refractory 



1. Thickness 

* 

* 

0.3 m 

2 . Thermal conductivity 

• 

1.3274 W/m°K 

3. Specific heat 


692.556 JAg.°K 

4 . Density 


3000 kg/m^ 

5. Emissivity 


0.9 

II. Solid 



A. Inlet temperature 

! 

366.66 °K 

B.Mass flow rate (dry) 

• 

2.996 kg/sec. 

C. Specific heat 

- • 

1122.54 JAg.°K 

D. Water in solid feed 

E. f Emissivity 

• 

18% (by weight) 
0.95 

III. Gas 



A. Outlet temperature 


9 35.77 °K 

B. Specific heat 


1214.7 JAg.°K 

C . Mass flow rate 


5.9075 kg/sec. 

D. Transmissivity 


0.5 

IV. Liauid-Water 



A. Latent heat of vaporisation 

: 

2.26xl0 6 JAg 

B. Specific heat (vapour) 

: 

2083.2 JAg-°K 
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V * Heat Transfer Coefficient 

A. Solid and refractory surface 

B. Outer wall and surrounding 
(only by convection) 


64 W/m 2 .°K 
10 W/tti 2 -K 


VI. Fill Angle 


: 90°. 
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Chapter 4 

Results and Discussion 


As has been mentioned earlier in the first 
chapter, the main objective of this work is the 
development of a kiln model which is more rigorous and 
improved than the simplified model of Sass [l}» Sass 
did not consider the circumferential variation of 
temperature in the refractory wall and his study is 
based mostly on various heat transfer correlations in 
contrast with the present work . 

The accuracy of the kiln simulation has been 
tested by comparing the present computer solution 
against that of Sass [l] and the actual kiln operating 
data. The set of data is obtained for a cement kiln 
from Sass [l ] and is the only set of complete kiln 
temperature-profile data available in literature. 
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The cement kiln temperature— profile data were 
obtained as part of a detailed analysis of cement kiln 
operation. The kiln was actually 185 ft. long, but 
since the present simulation neglects chemical reaction, 
only the 89.5 ft. portion of the kiln in which no chemical 
reaction occurs was used for comparison with the computer 
results . 

Figure 12 shows a comparison of temperature-profiles 
of the solid and the gas as a function of percent kiln 
length with those predicted by Sass [l] and with the 
actual measured data. Although some of the predicted 
temperatures are slightly low as compared to the other two, 
it is apparent that the simulation has accurately predicted 
the gas-temperature profile in the kiln, the fraction of 
the kiln used for water evaporation and the fraction used 
for heating the solids. 

However, the kiln length as predicted by the 
present model turns out to be 64 ft. as compared with 83.5 
ft. as predicted by Sass [ l] and 89.5 ft. which is the 
length of the actual kiln. 

The possible causes of deviation of the predicted 
results of the present work from those of Sass [l] and 
experimental data are the uncertainity in the input data to 
the computer program - for example. 



Temperature (°K) 


w 



FIG-12 AXIAL TEMPERATURE DISTRIBUTION 
OF SOLID AND GAS 
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(i) In Sass's paper no mention is made about the 
gas composition. The gas composition used in this work 
has been taken from a cement kiln hand book. Actually, 
the gas used in the kiln might be having different compo- 
sition . 

(ii) The fill-angle of 90° used here has been 
indirectly calculated from the data given by Sass. Sass 
has not explicitly indicated the actual fill-angle. 

(iii) Actual heat transfer coefficient from the 
inner wall of the kiln to the solid is not known. Only 
an estimate has been used. 

(iv) The actual ambient temperature outside the 
kiln is not stated in Sass’s paper. Also, the actual 
heat transfer coefficient from the outside of the kiln to 
the surroundings is not known. Here also, an estimate 
has been used. 

(v) No measurement of the solid— exit temperature 
or the gas-inlet temperatures are available. The values 
used are based on estimates made by operating personnel. 

Figure 13 shows the circumferential temperature 
variation of the inner kiln wall in the first axial 
segment of the kiln. The minimum temperature occurs at a 
point (© % 90°) where the wall emerges from the solid. 

After the wall comes out of the solid, it experiences a 

because it is then exposed to radiation 


rise in temperature 



Wall temperature ( K) 
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FIG-13 THE CIRCUMFERENTIAL TEMPERATURE 
DISTRIBUTION OF THE INNER KILN 
WALL IN THE FIRST AXIAL SEGMENT 

OF THE KILN 
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from the hot gas . The temperature rise continues until it 
is again covered by the solid (© = 360° or 0°) . However, 
it is clearly evident from the graph that the difference 
between the minimum temperature (742 °K) and the maximum 
temperature (748°K) is very small - which justifies the 
assumption by Sass that there is no circumferential 
temperature variation in the kiln wall. 

Figure 14 shows the radial temperature distribution 
of the kiln wall in the first axial segment of the kiln. 

It is clear that the temperature profile is almost linear. 



Wall temperature (°K) 
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/o length from inner wall 


FIG 14 ofthf A k D |! AL temperatur e distribution 
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Chapter 5 

Conclusion and Recommendation 

The present model developed for the non -reacting 
zone of a cement rotary kiln is an improvement over 
Sass 1 s model and can be used to design an actual rotary 
kiln used for drying and preheating of wet solids in 
many industries. Although the length of the kiln predicted 
by the present model shows a much lower value than that 
of the actual kiln, the axial gas and solid temperature 
distribution Vs. percent length of kiln matches well with 
the actual kiln profile. The main reason why the predicted 
kiln length doesn't match is the uncertainity in some of 
the input data. This can be overcome by checking the 
results of the theoretical model with a companion laboratory 
scale experimental set-up. For example, heat transfer 
coefficient from the wall to the solid, h , should be calculated 

C o 

from the data of an actual experimental work. Same can be said 
of the heat transfer coefficient from the outside kiln wall 
to the ambient. Lastly, future efforts should be made in 
the direction of the design of an entire cement kiln which 
includes- the reacting zone as well. 


Chapter 6 
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Appendix A 


A *° Calculatio n of Shane Factors 

The shape factors between the fifteen refractory 
surface elements and five solid surface elements are 
determined assuming the surfaces as infinite parallel 
strips . The equation for shape factors for infinite 
parallel strips was derived by Hottel [ 9 ] . Surfaces (1) 
and (2) in Fig. A.l are strips perpendicular to the XY 
plane and is measured from the normal to the surface ( 1 ) . 
The shape factor F ^_2 is given by 


1-2 


Sin 0 1 - sin 0" 
2 


(A.l) 


A detailed analysis of this shape factor calculation 
was done by Ghoshdastidar et al [ 4 ] . 


A * 1 Shape Factors between the Elements of the Kiln Wall 

Equation (A.l) is written in terms of the circum- 
ferential location of the two surfaces. For the notation 
shown in Fig. A. 2, the shape factor between surface (1) 
and ( 2 ) is given by 


„ _ A Q . Cosg 

F l -2 4 


(A. 2) 
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FIG-A-1 SHAPE FACTOR BETWEEN INFINITE 
PARELLEL. SURFACES 



FIG- A-2 SHAPE FACTOR NOMENCLATURE WHEN 
BOTH SURFACES ARE ON THE 
REFRACTORY WALL 
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where AO is the arc angle of surface (2) in radian. The 
angle B is defined as 


3 



(A. 3) 


where ©^ and ©^ are the circumferential position of 
surface elements (1) and (2) respectively with respect 
to the reference axis. For convenience, the reference 
axis is chosen perpendicular to the solid surface. 
Because of the absolute bracket in equation (A. 3), the 
shape factors F ^_2 an< ^ ^2-1 are e< l' ua l f° r the same AO 
as required by the reciprocity relations of shape 
factors . 


A. 2 Shape Factors between Elements of the Kiln Wall and 
Solid Surface 

In Fig. A. 3, surface (1) is a solid surface 
element and surface (2) is a refractory wall element. 
Surface (1) is located at a distance b from the reference 
axis and © 2 is the circumferential angle for surface 
(2) . Angle a is the fill angle of the solid material. 
Additional terms which are shown in Fig. A. 4 are needed 
to determine the shape factor. The shape factor Fj_ 2 is 
given by 

D Gos Cos A© 

F, , = 1 2 <*•«> 

1 4e 

where D is the inside diameter of the kiln»The angle 
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and f 2 are shown in Fig. A. 3 and are respectively 
the angles between the normals to surface ( 1 ) and 
(2) and the ray extending between the surfaces. The 
distance e is the length of this ray. Expressions 
for Y-jy 'F 2 and e in terms of parameters shown in 
Fig. A-4 were derived by Ghoshdastidar [4] and are as 
follows : 


^3 

= tan 

-1 (2b/D Cos 

a/2) 

(A. 5a) 

C 

= (D 

Cos a/2) /2 Cos ¥ 3 

(A. 5b) 

e 

CN 

0 

1 > 

II 

+ D 2 /4 -CD 

Cos (2m - © 2 + ^ 3)] 1/2 

(A. 5c) 

V 

2 

= Sin [ C Sin ( 2 it 

- © 2 + ^ 3 )/e 1 

(A. 5d) 

V 

i 

= e. 

- - TT 


(A. 5e) 
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FIG- A-3 SHAPE FACTORS NOMENCLATURE WHEN 
ONE SURFACE IS ON REFRACTORY AND 
THE OTHER IS ON SOLID 



FIG-A-4 PARAMETERS FOR SHAPE FACTOR 
EQUATION 
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Appendix B 

Calculation of the Emissivity of the Gas 


For the calculation of the emissivity of the hot 
gas, we require the composition of the gas. Of the gases 
encountered in heat transfer equipments, carbon monoxide, 
the hydrocarbons, water vapour, carbon dioxide, sulpher 
dioxide, ammonia, hydrogen chloride and the alcohols 
are among those with emission bands of sufficient magnitude 
to merit consideration. Gases with symmetrical molecules, 
hydrogen, oxygen, nitrogen etc. have been found not to 
show absorption bands in those wavelength regions of 
importance in radiant heat transmission at temperatures 
met in industrial practice. 

The composition (by volume) of hot gas for a 
typical cement rotary kiln was obtained from Ing. Paul 
Weber [5], and is as follows: 

C0 2 - 15.897* 

0 2 - 1.418% 

N - 41 . 784%'_ 

2 

H 2 0 - 40 .9% 2 

Thus we see that only C0 2 and water vapour are 
the gases having any effect on emissivity. For calculation 
of emissivity, we are using the method suggested by Holman 
[6 ], using mean beam length, Le. 
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Since mean beam length information for this particular 
geometry is not available, as per Holman [e], we can use 
the following approximation 

Le = 3 - 6 a (B.l) 


where V is the total volume of the gas and A is the total 
surface area. For the fill angle of a = 90° in the 
present work, we can see that 


Le = 3.6 


X"R 2 

4 


1 r 2 
2 R 


3 

2 * R 


+ /2 


i.e. Le — 1 .6783 R. 

where R is the internal radius. 

If the pressure inside the kiln is assumed to 
be 1 atmosphere, then partial pressure of C0 2 and water 
vapour can be obtained as 0.1589 atm. and 0.409 atm. 
respectively. Radius of the kiln is 3.25 feet. 

Using the above values and the method given in 
[ 6 ] and converting the temperature scale into degree 
Kelvin, we can tabulate the emissivities as in Table 2. 
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Table 2 . Emisslvitv of Hot Gas at Different Temperature 
Ranges 


Temperature 

range °K 

Emissivity 

From 

To 

273.0 

347.0 

0.6304 

347.0 

486.0 

0.5524 

486 .0 

625.0 

0.5358 

625.0 

764.0 

0.5408 

764.0 

903.0 

0.5210 

903.0 

1042.0 

0.5028 

1042.0 

1181.0 

0.4680 

1181.0 

1320.0 

0.4170 

1320.0 

1459.0 

0.3838 

1459.0 

1598.0 

0.3506 

1598.0 

1737.0 

0.3202 
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Appendix C 

C .0 Finite Difference Scheme 


A. Inside Surface Exposed to Radiation from the 
Hot Gas 

Refer Fig. C-l. Here q is the radiation heat 
input into the element on the refractory walls from 
the hot gas . As per the coordinate system already 
explained, we can see that 


AY = 


. 2tt _ a 
15 


where is the fill angle in radian 


_ , ThicknesS ' _ 0.3 

AX - l 1Q ) 


10 


= 0.03 m 


Then finite difference equation can be written as 

, J >+1 
AX 
2 


-K AY. AZ 


rpP _ipP 

m.,n m+l,n 


AX 


+ q = p C AY AZ 


mP mP 

+ p s cM 4Y4Z . 


T'" - 

m , n m , n 

At 


(C.l) 


B. Inside Surface Covered by Solid 


Refer Fig. C.2. Here, the area covered by the solid 
is divided into five elements. So 





47 


AY = R . f 
AX = 0 .03 tti. 


In this case, the finite difference equation can 
be written as follows: 


T 9 - T 9 
m.n m+1 . n 

AX 


-K.AYAZ . = p C -r^AY.AZ 


AX. 
s - 2 


T P+1 - T P 
m,n rn,n 

At 


+ h e (T p -T ) AY.AZ + p C 
cs m,n s K s 


AX 


.AY.AZ. U 


T 9 - T 9 - 
m, n m, n-1 

AY 


(C-2) 


C. Interior Point 

Refer Fig. C-3, For a general interior point. 
Ay can be calculated as 


AY = [r + (K— 1) AX J 


2ir - a 
. 15 


where K is the grid number in X direction and 


AX = 


Thickness 

10 


= 0.03 m. 


Then finite difference equation can be written as 


follows : 
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TP + T P - 2T P 1 tP +1 _ 1^ 

w+l i n m-1 , n *n, n — 1 m.n m,n 

2 a At 

(AX) rf 


+ 



ijP _ i<P 

m, p m.n-1 

AY 


(C. 3) 


where a 

-rf 

material . 


is the thermal diffusivity of the refractory 


D. Outer Surface Exposed to Ambient Air 

Refer Fig. C.4. If hr is the combined radiation 
and convection heat transfer coefficient for the outer 
surface, then the finite difference equation can be 
written as follows: 


T p+1 _ tP 

„ . v ,„ x ^m-l,n _ ^i,n _ r AX „ ?i f n Hle.P. 

-K.AYAZ ^ “ P s C * 2 * AY * AZ Ax 


(T? « - t£ 


AV tP - T 5 

+ h (tP - t )ay.az +p c. .ay.az. u. ^^ -' 7 m ' p - L 

_r inf s ^ AY 


where T is the ambient temperature in degree Kelvin. 

In the above cases, the circumferential velocity, 

ttDS 

U, is calculated using the formula, U = — q- where S 
is the number of revolution per minute of the kiln. 


C.l Method of Solution 

Equations (C.l), (C.2), (C.3) and (C.4) are 

solved simultaneously at each time step At using explicit 
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method. A time step of At = 1 second was taken. Since 
the initial temperature distribution of the refractory 
wall was not known, an arbitrary temperature was assumed 
at all grid points in the refractory wall. The solution 
converges when there is no further change in the temper- 
atures at each grid point as x -*■ » . This temperature 
distribution represents the steady state temperature 
distribution of the wall. 
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CENTRAL LIBRARY 

1 i.T.. KANPUR 


Ace. No. AXQ&QM 
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*************** + ***********#******* + ************♦**¥***********■$:**** 

program for simulation cjf near transfer in ccmf.nt rutark ku.m 

********************************* *********************************** 

DTMRNSTQN FC?0*20) , T'U20l r TTH( 70 3 # 1( 5) , f-S( ?t> 3 

DIMENSION OR 120) , AC20) , FG (7o) ,Te'V( 20) , jEMPTf 2?,H ) -OCC5) 

dimension tswi700),tem£gC2003 

0PFN(UNIT=2S,DFvTCF='DSK' ,F] r,£='r.)R2S.D ft r ' ) 
npF^f ljNlT = ?6 / DFVTCE='DSK',FH,E='rOR26.DR I*' ) 
****♦************************+************♦*****♦****************« 

definition of terms used 

****************************************************************** 

D - - ■ ~ Inner diameter of the Kiln ( M ) 

Alpha - - mi angle of the solid (Radians) 

T FMpG - - Temperature of the gas Chenpee lO 

DpTi r - - Time interval In seconds 

SiGma - - Stefen Boitzman constant 

R 0 - - - Density of retractorv KG/C U . M 

CpR - - - Specific hent of, rofrarLn r v (j/kG,k> 

cpS - - - .Specific neat of soli j CJ/KG-K) 

CpG - - - Specific heat of gas td/KG.iO 

Hc-S - - - Convective heat transfer coef f C jnsl de) (N/Sq m.k) 

Hrs - - - Combined heat transfer coef t »( Outside) C w /Sq A * K > 

Ts w - » - Temperature of solid (ocorep k 1 
TjNf - - Ambient temperature fnegreo K) 

PK" - - - Conductivity n f refractory (W/M.K) 

ThICF - - Tnic'cness ot refractorv w-ii i (M) 
u - - - - Ctrcupf erent tal velocity of solid CM/Sec) 
tqwg - - Transmissivity of oas 
FpSG - - Fmtssivity of a as 

SmASS - - Mass flow rate of drv solid (K'>/Rec) 

GMASS - - Masslw rate of hot. gas IKG/Sec) 

CM A x - - Ouantlfy of water tobe evangrafed (KG/SeC) 
it i T , ft T - ~ Latent neat, of vaporisation of water (l/^Gi 
HtjDsR - - Total heat loss from onLe r solace (N) 

Or - - - Radiation heat transfer Cwi 
0C - " - Convection heat transfer (N) 

Tempt - - Temperature of wall f Degree Kj 
F(T,,T)- - Shape factor 

: * * * * * / 

h = 1 .Oij 

T E M pG ( t )=9lS.77 
D Z = 1 * U 

T 3 W ( 1 )=3bh.bA 
T5 h AX=077.1 



r. 

r 

r 


FnAX=0.5l93 
r> J.riMArS . 6b9£-8 

Rn= 30 oo.o 

CpR = 692.5;5r> 
nc=!0.0 
H C -R = 6 4 . o 
TllTCK= f >.3 
”=~u.iRun 
nei. f=i . 

T 'J*'S=9.5 
H£T.ATs2.26E*6 
CpV=?0«3.2 
OpS = t l?2.54 
p r*3. 14 15926 
ALR;iA = Pl/2. 

PHjT.TMsf 2.*Pi-AfiP»lA)/J 5. 

p =n/ 2 . 

**********************»***♦•# + * + *** 
r AT-CMLjATTuM OK SMAPK F'ACTup 

R C 1 )=-R.tR + STiJ(AF»F’'»A/2. )/]n, 
RC2) S H(1)/?. 

P>c ^ * 0 - 

p(4J=-t ,*Bf 2) 
nc r O=-t .*n(l) 
r>n 666 x=i,2n 


)■ F t * t f 


TO 1 1 1 
)*AI,PMA/1o. 


Ill 

ll? 

666 


F, = T 

TF(i;.G^.6) G f1 
Tiur) = f2.*r.-t 
Go jpn U2 
T l'f(ll=ALPHA|-r2-^f 
riVJ rTN’iE 

TrH(r) = rH(T j-a».PM’A/2 . 
rri^TTM^E 

r>U 777 1 = 1,20 
>■>0 fifty J=1,20 
f f ( r.GE.6) on 
Tp (I. OF.. 6) oy 


t *r*Kf, I 0/7. 


F( r, jj=o.o 


T n 

TO 


im 

1H2 


1 1 1 2 


f cg' t 5 = o - 
On yn Ry» 

R:;T3 = ArA'H2.*HCi)/N)*rn- 1 fAr 1 nnA//. ) ) ) 

r = nfC06( Al.PHA/2. l/f 2, *0;). r >( o S T j 1 ) 

p:=S 3 R i' ( c*c+o*i>/ 4 .-r’ + o»roo ( 2. krt-ti m { .n *.ps i 3 > ■> 

n 5 t 2 = as I fj c r * 6 r m ( 7 . * p t-ttm ( , 7 ) +ps r i j / e 1 

p STlsfTHCj)-P.';l2-Pt 

Ff T ,0 ) =u*cnsf PS J 1 ) #cn,sf po 1 7) I'OC’IjT-I / t 1 . »■(.; ) 

F(0, T )=u.4*STiJf hF'POA/2. )H ( l , J 1 / r>u>. V N 
0;j rn ftyrt 

Tp( J.I.F.5) on T n ! t Rn 


1 1 1 1 
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bet a= A f\s ( (absctth(T)-i t h f o) )-pt ) / i . ) 

F(T, J)=DELTH*C0s(HF:rA)/4. 
flg8 CGNTTNUe 

ill continue 

c 

c verification of shape factors 

c - - — — - - ... ... 

oo 2222 1 = 1 ,20 

FSCI)=0.0 

no 2221 J=i,2o 

FS(I)=FS(I) + F(T,,I) 

2221 CoHrlN'lE 

2222 ^ COMITNUE 

c t************************** ********** **+********* 

C WRT rB( 25.999 ) C(F(T,d) ,T=1 ,30) ,J=l ,20) 

Cylo F0R«MC(20fF6.3))/) 

C WRTrE(25,998)(Fs(I),I=l,20) 

C998 FqRmAiC () 0X,F8.5>/> 

00 21 31=1 ,20 
TFCr.I^.S) GOTO 2H 
A (T)=R*D7,*DElrrH 

Gqt 3 2 1 3 

2 1 4 A(T)=0.4^R*L)7,tSl''FCAL D h A/2.1 

2 i 3 continue 

AG = 6. l266*R*nz 
T I N F- 300 • 

0K = U 3274 
ox=thick/io. 

OQ 215 1=1,20 
2 1 5 F G CI)=ACI)/AG 

C ----- — — — — - 

c initial assumption uf Mali, gpio temp. Cuporee kj 

' "MEADC47"696)f C T^MpT ( l",T) , J=7, 77) , i=1 ,2 n ) 

696 FuRMMC(l5X,1lC2X,F0.Sl)/) 

fm=o.o 

Hr,oss=o.o 

no 100 m=i , 30 

“Hi TE ( 2b . 992 ) ( M TSW ( M ) ,tFMPC.( O ) 

992 FuPMATClOx,l3,sj(,FO. : »/»X,FM..3) 

00 300 N=1 ,3000 

C ------------------ 

c CaT-COLATToM of w A l>l. TpMp. U5T-JC FJ m |TF: orFFfc-PpMCF MF.Tt.ni) 

c 

TDEljSO.O 
00 301 J=l,11 

on 302 1=1,20 

IF (.J.NE.l) GOTO 375 

00 2H K=l,5 

2 1 1 T5MpCKl=TsW(M) 

00 212 K = 6 , 2 0 
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51 ? 
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521 
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tf 
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GuTf) 
nnTo 
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GfjTd 

GOTO 

corn 
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GOTO 


5 1 0 
SI 1 
SI 2 
SI 3 
ST 4 
SJS 
sir, 
517 
Sin 
s i 9 
5 20 


(TE«PG(M) ,LE. 374 . ) 

(T£MPGC 0 ) ,l,F. 4 bR. ) 

(TEMpG(M) .f.E. 625 . ) 

(TE 5 PGCH).|,F. 7 £» 4 .) 

(TEMpfi(H) .LF, 9 (; 3 . ) 

(TEMPG(M) .I,E. 1042 .) 

(TE^pG(M) .i.R. 1 181 .) 

(TFMpG(M) ,i,F. 1 120. ) 

(TE*lpG(M) ,UK. 1 45 ‘j.) 

(TE'lpG(M) . i,F. 1 58 H. ) 

(TEMpG( 5 ).GT. 150 f».) 

FPSG = 0.6 3 C '4 

Go T n 52 ] 

Fpr>r, = 0 . 5 S 24 
goto 521 
F PSG^O.STSS 

fqto 521 
fps:,=o. 54 ofi 
GjT 3 521 
EPS 5 * 0 . 5 ? 10 
Gf)T 3 5 ? 1 
FPS 3 = 0 . 502 « 

Goto 521 
Fp.53so.46uo 
goto 521 

EPS8=0.4l 70 
GU'Tf) 571 
EPSK=0, JO )« 

gutjj 521 

EPSS=0.3536 

goto 571 

FP8 3 = 0.37 07 
FfjM 1 ’TN'lF 

PpcftGtFGC I ) *EPS(;»STG’lft*( T t ;MpGiM )* v 4 ) 

AKKsO.O 
n !J 717 ?j=l , 2 f > 

A KF=A KK + A ( U) *F C U , .O * ToWf;t 3 1 G ** i\+ ( TKM p f L ) * * 4 ) 

C 1 )H rT.^'JE 

FtjA I = A Cl ) *K 1 GH A* fTFMPl T ) * + 4 ) 

Of< f I ) =Pp + A«K-EOAT 

C ******* *+»** + *# + . t******** 4 .**: + :$*:};.M***** 4 ’* 4 *f*ic*** 

375 G !l M r T M 0 E 

C WRI rF.(?5,907) ORU ) , T =1 ,?U) 

C9Q7 r fJ RMArC( 7 oX,Fl 5 . 4 )/) 

TEMPTI?! , J)=TEMpTt 1 ..I) 

r>Y=(R+- ( J"*l ) *DX) ♦OKIiT'O 
TMFrfSO . 0 
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7. ( M*Dr!jT3 /OY 

Z2 = c 2 . *pK»0li;r. , r) / C R0*CpF?*f)X*l)X } 
rpc j.fjKj ) cnxn 37 fi 
TFd.KO.I) goto 377 

TP ( T . r . E . 5 ) GOTO 37 n 


gktd 0 0 x m t g nr; the iumfr waij, 


Z i= c 2 . * dhlt J / ( RO*CPR * l')K * r>y *i)Z ) 

TT2sTE^pTCT.1)*(1. +7.1-72) 

T.jFrfsrT2 + Z2*TEMpT(T r ?)-Zl1tTEMPT(T-1 ,1)3 Zi+ORf [)*n.9ij 
flilTO 390 

7 l s ( 7 . *DP:lT + MC>3 ) / ( Rf)*C°R*OX) 

7l'3=T E wpTCT,1 )*(t.+Zl-7.2-Z4) 

Ti]F<tf=rT3 + Z7+T E MpT ( T , 7 ) -7 1 *TF«PT C T - 1 , 1 ) + ?,4tT5 ! ; f M ) 

CrjT3 3-9o 

Z4=( ?.*DEl,T + HCS) / lRO*Cd**OX) 

TT4rT E MpT ( 1 , 1 ) + ( 1 .+71-7.2-7'!) 

T[]p^ = xT4 + ZZ + Tf’ v pT(, 1 ? ) - 71 *T E M pT ( 70 , 1 1 +7 4 + '|'S w C n ) 

G n T ;j 3 9 0 

T F fJ.'RQ.H) GOTO 3/9 
TF fT.FQ.l) G»To 390 


INTERIOR GPIO FUTflTS 


Tl’r ( Tt; M pT ( I ,1+1 ) 4 TPMf’TC I , J - l ) ) + 7.7/7. 

T0F^ = IFMPT( 1 , J ) + Cl . +7, 1 -77 ) 4 r T-7. 1 * i'P-M M r < I - 1 , J ) 
OOTD 390 

Tr5sflFMPT( 1 , J + l )+TR*1pT(l J-J )) *7.2/2. 
TfiPWsTEMPTf 1 , J)*f 1 . + 7 1 -7,2) + TTS-fcl * TfiUpT f 2y , J) 
goto 3 q o 


GhTo POINTS ntl THE 11 IITFK F A | , f . 


MCI =(TFHPTf 7. , J) + TI»F)*f TFHrrC I , J) **?.* HNF**2) 

HRSsHC+HCl*STGMA*U.9 

7.[i~ (?.* OKI.' T* MRS) /(8(lt C P R*L>X ) 

IF (T.FQ.l) G0 T 0 3«1 
Tt 6 = T E » PTcT , 1 1 ) * ( ) . + 7 . | - 7 , 2 ~ 7 .b ) 

T I ( P W = T T Q + 7 . 2 * T f j M P T ( 1 , 1.0 ) -7,) *1 FMP'I f 1 - 1 , j 1 ) + Z5*T J»F 
G()TJ 39 0 

Ty 7 =T E M pT ( 1 \ 1)*M . + 7.J-JS7-7.5) 

Ti!F^=rT7 + z?*T|.-:»pT( 1 , 1 i.) ) - 7. 1 * T F M P rf 20, 1 )) + 7C>*TTbF 

run n nme 

T|)F'Li = TOFF yAR.S ( T fJFU-TFt'D'r f j , J ) ) 

T|VMpT( r iJ )= TNE w 
Cn fI FTN n K 
r;inrT,\UE 

TFf TDEI 4 .LF. . 1 ) GHTO 7u9 
rjMTlNOE 

continue 



r 


096 

c 

r 

c ** 


4 U 9 


r. 

c 

c 


410 


} l* 

C - 

C 

c - 


412 

C - 

c 

c - 


C 

c - 


4 30 


r 

r: 

c 


r. - 
r. 

c ~ 

460 


TF (M.GE*3) GllTO 409 
Wim’F(25, 9 q 6) C CTRMPTCI, J) , d = t 
FQR4ATC C'5X,11 t2X # Pa,31)/) 


CAliCHL^HOM HF HFAT MRS 


n,T = i , 2 0) 


Hbnss=HLnss+HRS»f)2*PT)»r)* ( vfmp r ( 1 , 1 1 j-j'TnF) 


njNvPCTjnM rn sor.in from refractory w ^ r, fj 


no 410 1 = 1,5 

0Cri)=HCS¥A(6)*(Tfc;MpT(T, 1 }-TR;V f M) ) 

SlGyCsO-.O 

00 HI 1 = 1,5 

SiGjC=SlGor+OCCl> 

C()M f XNtJE 


R A n t 4 1 TQM TO GOT, ID FROM G^S 


SlGiRsO.o 

00 412 1 = 1,5 

Gx 0 3 R = RlGoR + O .97 5 * 0 Rm 

gunitnoe 

SI.G 3 RT=STCJ 0 R 


ToTAO RAOJ.ATTOM UFA I.’ TRa’U-jFKR PRO i 3 AS 


Of) 413 I =(,,20 

Sl.noRT = STOORT + (>. f>5*QP ( T ) 

CrjfJTTNME 


r Al'C r, r,ATlOO op SOl.TO AM|) (JAR I’Frp. 1 M T q K Hp.XT RF(J Op rI 1 ’.. 


Tp (FM.GF.FMAX) Gu t O 4 5f, 
tf (TrW(D) .po. 377 . •» goto -Ho 

TsW(M + t ) = T G w f fO ) rSTGOR + sTGOCl/f (R > A 3 .3 *- F. i A X ) * 0 P 5 ) 
TF (TSW(M *1 J.bT. 373 .) G()T tj 4 -Ip' 

T 3 ’■( ( M + 1 ) = 37i.O 

goto 440 

I F (PM.MP.O.-O) GOTO InO 


Pi, J - r.pr'(,TH of preheat j/'g x.nr.jF 


M 1 = V) - 1 

Poi =mi*ut: 


OfTP - AMOijM-p 1 jF OATpP Fv-APMOATFI) j. m t .{ T s oFiO'.pMT 
OpO = C S T GOp + 3 T GQC 1 / 0 1 0 at" 

R T = ( G * l A S G + F M A X - F . '* ) / C G >1 ASR) FM A X - F'f -O i" ) 



7 


4 40 
450 

C - 

c 

c - 

4 ft 0 
4gO 

c - 

c. 

c - 


c 

c. 

c 


99 


492 

loo 

5 y 0 

995 

991 


a m t r=r> t nKM*n*v*('rnMPR f 0 ) - 174 . 0 ) 

Tl^'pCJc'l + l ) = T F cl P d f Ml ♦ 0 T +• A M I’T/ ( cfMAS,9 + FriAx-|'f' j trPf 

TsW(M+t )=37.U0 

r.QT'j 9 9 

Tf.;f'pG(fN 1 )=TFMPnf rO+K.lRo" !•/ 1 ( |.;m\X ) t(’PC) 

f’tlT 9 9 ° 

t p’ tTgwiH) . i;Q.rn. 1 nnrn 4^0 

ROT a 490 


P(. t ? - IjF.MRTH UP’ RVApnRAl'TNU 7,'.|Pp; 


PIj 2 s(m 2 -M 1 )tnz 

T 5 W ( M + 1 ) =T RW ( M ) + C STROk + r> l «0 C V f S"’ A Sii * CPS J 
T f M p G ( M 1 1 ) = T F. M P R ( M 1 + 9 J R q ft | / ( C, ’,] A 'J 9 -I F P R ) 


PO 3 ” LENGTH OF PORTMEAtTNR ZUPpJ 


N i - M -m 2 
P| ( *sM 3 *DZ 


pi. - r n t a i , i.FNRiu of riiE kiln 

p ^ ij[ lj ^ 

TP~(T5W( 1 ) .GF. TS'-IAT ) G n m 5oO 
CON flNHfc; 

PQ 492 1 — 1 . r 2 0 
Hfj 492 J = 1,U 
Te:MPT( r,U)=TFHPTfI , JH9S. 

r,ri T'T'iof: 

w R T l’F(26»9ft21 ( T 5 W ( ^ i 1 ) , TFmPR f 4 + 1 ) ) 

fo^a c c i*x, f<>. i, sx , f»>. ) > 

^RTfF(20,9P5) CPM , P»-2 , »>»• i , PI. ) 
FC)9MAX( ( I OX. 1 FI 2. 3)/) 

PRT rF(2b,99i) HunsiR 
FQRv|At(20X,F1 5.3) 

Si' OP 
Ffjn 



